Body state estimation of a vehicle

ABSTRACT

The present invention features a system and method for estimating body states of a vehicle. The system includes at least two sensors mounted to the vehicle. The sensors generate measured vehicle state signals corresponding to the dynamics of the vehicle. A signal adjuster transforms the measured vehicle states from a sensor coordinate system to a body coordinate system associated with the vehicle. A filter receives the transformed measured vehicle states from the signal adjuster and processes the measured signals into state estimates of the vehicle, such as, for example, the lateral velocity, yaw rate, roll angle, and roll rate of the vehicle.

BACKGROUND

This invention relates to a system and method of estimating body states of a vehicle.

Dynamic control systems have been recently introduced in automotive vehicles for measuring the body states of the vehicle and controlling the dynamics of the vehicle based on the measured body states. For example, certain dynamic stability control systems know broadly as control systems compare the desired direction of the vehicle based on the steering wheel angle, the direction of travel and other inputs, and control the yaw of the vehicle by controlling the braking effort at the various wheels of the vehicle. By regulating the amount of braking torque applied to each wheel, the desired direction of travel may be maintained. Commercial examples of such systems are known as dynamic stability program (DSP) or electronic stability program (ESP) systems.

Other systems measure vehicle characteristics to prevent vehicle rollover and for tilt control (or body roll). Tilt control maintains the vehicle body on a plane or nearly on a plane parallel to the road surface, and rollover control maintains the vehicle wheels on the road surface. Certain systems use a combination of yaw control and tilt control to maintain the vehicle body horizontal while turning. Commercial examples of these systems are known as active rollover prevention (ARP) and rollover stability control (RSC) systems.

Typically, such control systems referred here collectively as dynamic stability control systems use dedicated sensors that measure the yaw or roll of the vehicle. However, yaw rate and roll rate sensors are costly. Therefore, it would be desirable to use a general sensor to measure any body state of the vehicles, that is, a sensor that is not necessarily dedicated to measuring the roll or yaw of the vehicle.

BRIEF SUMMARY OF THE INVENTION

In general, the present invention features a system and method for estimating body states of a vehicle. The system includes at least two sensors mounted to the vehicle. The sensors generate measured signals corresponding to the dynamic state of the vehicle. A signal adjuster or signal conditioner transforms the measured vehicle states from a sensor coordinate system to a body coordinate system associated with the vehicle. A filter receives the transformed measured vehicle states from the signal adjuster and processes the measured signals into state estimates of the vehicle, such as, for example, the lateral velocity, yaw rate, roll angle, and roll rate of the vehicle.

The filter may include a model of the vehicle dynamics and a model of the sensors such that the states estimates are based on the transformed measured signals and the models of the vehicle dynamics and sensors. The filter may also include an estimator implemented with an algorithm that processes the transformed measured vehicle states and the models of the vehicle dynamics and sensors and generates the state estimates.

The present invention enables measuring the body states of a vehicle with various types of sensors that may not be as costly as dedicated roll or yaw rate sensors. For example, the sensors may all be linear accelerometers. However, in some implementations, it may be desirable to use an angular rate sensor in combination with linear accelerometers.

Other features and advantages will be apparent from the following drawings, detailed description and claims.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 depicts a block diagram of the processing of the vehicle states in accordance with the invention.

FIG. 2 depicts a general array of sensors for measuring body states of a vehicle.

DETAILED DESCRIPTION

In accordance with an embodiment of the invention, FIG. 1 illustrates a system 10 that measures the vehicle states of a vehicle identified as block 12. Specifically, the system 10 includes a plurality of sensors 14 that measure signals which contain parts related to components of the vehicle states of the vehicle dynamics 16 produced, for example, when the angle of the steering wheel δ is changed.

The system 10 also includes a signal conditioner or adjuster 18 that receives measured signals from the sensors 14 and a filter 20 that receives the adjusted signals from the signal adjuster 18. In certain embodiments, the filter 20 is a Kalman filter including a model of the vehicle dynamics 22 and a model of the sensors 24. These models are described below in greater detail.

The signal adjuster 18 and the sensor model 24, which incorporates the model of the vehicle dynamics 22, provide inputs to an estimator 26. An algorithm with a feed back loop 28 is implemented in the estimator 26 to process the transformed signals with the models of the vehicle dynamics and the sensors. The output from the estimator 26 is the state estimates {right arrow over (x)}_(v). The body states estimates may include the roll angle, roll rate, yaw rate, and lateral velocity, as well as other body states.

In some embodiments, the sensors 14 measure the linear acceleration at a particular location where the sensor is mounted to the vehicle. When the sensors are not aligned in a plane perpendicular to the axis of interest, the measured values contain biases proportional to the angular rates about other axes. Similarly, when the measurement axes of the sensing devices are not coincident, the measured values contain biases proportional to the angular acceleration about other axes. Moreover, when the measurement axes of the sensing devices are not coincident and are not mounted along a body reference axis, the measured values contain unique gravity biases dependent upon the difference in mounting angle of the sensors and the body lean angle of the vehicle.

To address these biases, a general implementation of the system 10 can be employed as illustrated in FIG. 2. Here the sensors 14 (identified individually as S₁ and S₂) are in known and fixed positions on the vehicle body 12 and the orientation of the measurement axes of the sensors S₁ and S₂ are known and fixed. Specifically, the location and orientation of a sensor S_(i) is provided by the relation P_(i)(x_(i), y_(i), z_(i), θ_(i), χ_(i), φ_(i)),   (1) where x_(i), y_(i), z_(i) are the space coordinates of the sensor S_(i), θ_(i) is the sensor yaw angle, that is, the orientation of the sensor's measurement axis in the X_(B), Y_(B) plane with respect to the X_(B) axis, χ_(i) is the sensor pitch angle, that is, the orientation of the sensor's measurement axis with respect to the X_(B), Y_(B) plane, and φ_(i) is the sensor roll angle, which is the rotation about the respective measurement axis.

The sensors S_(i) measure the linear acceleration at the location Pi, namely, {right arrow over (α)}_(i)={right arrow over (m)}_(i)·|m_(i)|=[α_(xi), α_(yi), α_(zi)]^(T), where {right arrow over (m)}_(i) is the unit vector along the measurement axis, and |m_(i)| is the magnitude of the acceleration along the measurement axis.

Since the acceleration {right arrow over (α)}_(i) measured by the sensor S_(i) is the acceleration in the sensor coordinate system, the measured accelerations are transferred to a body coordinate system. In certain embodiments, it is assumed that in an array of single axis accelerometers each accelerometer has a measurement axis referred to as the x_(sensor) axis. Accordingly, the transformation from the sensor coordinate system to the body coordinate system is provided by the expression $\begin{matrix} {\begin{matrix} {{{\overset{\rightharpoonup}{a}}_{i} \times {\overset{\_}{Body}}_{i}} = {{{\overset{\rightharpoonup}{a}}_{i}\begin{bmatrix} x_{{body},i} \\ y_{{body},i} \\ z_{{body},i} \end{bmatrix}} = \begin{bmatrix} a_{x,{body}} \\ a_{y,{body}} \\ a_{z,{body}} \end{bmatrix}}} \\ {{{where}\quad{\overset{\_}{Body}}_{i}} = \begin{bmatrix} x_{{body},i} \\ y_{{body},i} \\ z_{{body},i} \end{bmatrix}} \\ {= {\begin{bmatrix} {\theta_{i}^{c}\quad\chi_{i}^{c}} & {{{- \theta_{i}^{s}}\phi_{i}^{c}} - {\theta_{i}^{c}\quad\chi_{i}^{s}\phi_{i}^{s}}} & {{\theta_{i}^{s}\phi_{i}^{s}} + {\theta_{i}^{c}\quad\chi_{i}^{s}\phi_{i}^{c}}} \\ {\theta_{i}^{s}\quad\chi_{i}^{c}} & {{\theta_{i}^{c}\phi_{i}^{c}} + {\theta_{i}^{s}\quad\chi_{i}^{s}\phi_{i}^{s}}} & {{{- \theta_{i}^{c}}\phi_{i}^{s}} - {\theta_{i}^{s}\quad\chi_{i}^{s}\phi_{i}^{c}}} \\ {\quad\chi_{i}^{s}} & {\quad{\chi_{i}^{c}\phi_{i}^{s}}} & {\quad{\chi_{i}^{c}\phi_{i}^{c}}} \end{bmatrix} \cdot}} \\ {\begin{bmatrix} x_{sensor} \\ y_{sensor} \\ z_{sensor} \end{bmatrix}} \end{matrix}{where}{\,_{\_ c}{= {\cos(\_)}}}{\,_{\_ s}{= {\sin\quad(\_)}}}{\theta_{i} = {{sensor\_ yaw}{\_ angle}}}{\chi_{i} = {{sensor\_ pitch}{\_ angle}}}{\phi_{i} = {{sensor\_ roll}{\_ angle}}}} & (2) \end{matrix}$ and [x_(sensor) y_(sensor) z_(sensor)]^(T)=[1 0 0]^(T), since x_(sensor) is assumed to be the measurement axis for each of the single axis accelerometers.

Note that the transformation identified in Equation (2) is typically performed in the signal adjuster 18 (FIG. 1). The signal adjuster 18 may also provide a DC bias offset compensation to compensate for the biases discussed above.

Regarding the Kalman Filter 20, the model of the vehicle dynamics 22 for a state vector {right arrow over (x)}_(v) =[{dot over (y)} _(v) r _(v)θ_(v){dot over (θ)}_(v)]^(T)   (3) is provided by the expression $\begin{matrix} {{{\overset{.}{\overset{->}{x}}}_{v} = {{A \cdot {\overset{\rightharpoonup}{x}}_{v}} + {B \cdot \overset{\rightharpoonup}{u}}}}\quad} & (4) \\ {{where}\quad{\quad{\begin{bmatrix} {\overset{¨}{y}}_{v} \\ {\overset{.}{r}}_{v} \\ {\overset{.}{\theta}}_{v} \\ {\overset{¨}{\theta}}_{v} \end{bmatrix} = {\begin{bmatrix} {- \frac{C_{F} + C_{R}}{mu}} & {\frac{{C_{R}b} - {C_{F}a}}{mu} - u} & 0 & 0 \\ \frac{{C_{R}b} - {C_{F}a}}{I_{z}u} & \frac{{{- C_{F}}a^{2}} + {C_{R}b^{2}}}{I_{z}u} & 0 & 0 \\ 0 & 0 & 0 & 1 \\ {- \frac{h}{I_{x}u}} & \frac{h\left( {{C_{R}b} - {C_{F}a} - {mu}^{2}} \right)}{I_{x}} & {- \frac{K}{I_{x}}} & {- \frac{C}{I_{x}}} \end{bmatrix}\quad{\quad{{\begin{bmatrix} {\overset{.}{y}}_{v} \\ r_{v} \\ \theta_{v} \\ {\overset{.}{\theta}}_{v} \end{bmatrix} + {{\begin{bmatrix} \frac{C_{F}}{m} & 0 \\ \frac{C_{F}a}{I_{z}} & 0 \\ 0 & 0 \\ \frac{C_{F}}{m} & 0 \end{bmatrix}\quad\begin{bmatrix} \delta \\ g \end{bmatrix}}\quad{and}\quad{where}{\overset{.}{y}}_{v}}} = {{{lateral}\quad{velocity}\quad{of}\quad{the}\quad{vehicle}r} = {{{yaw}\quad{rate}\quad{of}\quad{the}\quad{vehicle}\theta_{v}} = {{{roll}\quad{angle}\quad{of}\quad{the}\quad{vehicle}{\overset{.}{\theta}}_{v}} = {{{roll}\quad{rate}\quad{of}\quad{the}\quad{vehicle}C_{F}} = {{{cornering}\quad{stiffness}\quad{of}\quad{the}\quad{front}\quad{axle}C_{R}} = {{{cornering}\quad{stiffness}\quad{of}\quad{the}\quad{rear}\quad{axle}a} = {{{distance}\quad{from}\quad{center}\quad{of}\quad{gravity}\quad{to}\quad{the}\quad{front}\quad{axle}b} = {{{distance}\quad{from}\quad{center}\quad{of}\quad{gravity}\quad{to}\quad{the}\quad{rear}\quad{axle}m} = {{{mass}\quad{of}\quad{the}\quad{vehicle}h} = {{{height}\quad{of}\quad{the}\quad{center}\quad{of}\quad{gravity}\quad{above}\quad{the}\quad{roll}\quad{axis}\quad\text{}I_{z}} = {{{yaw}\quad{moment}\quad{of}\quad{inertia}I_{x}} = {{{roll}\quad{moment}\quad{of}\quad{inertia}C} = {{{vehicle}\quad{roll}\quad{dampening}K} = {{{vehicle}\quad{roll}\quad{stiffness}u} = {{{longitudinal}\quad{vehicle}\quad{speed}\delta} = {{{steering}\quad{angle}\quad{of}\quad{the}\quad{tires}g} = {{{gravitational}\quad{acceleration}\overset{.}{*}} = {{\frac{\mathbb{d}}{\mathbb{d}t}*{and}\quad\overset{..}{*}} = \frac{\mathbb{d}^{2}}{\mathbb{d}t^{2}}}}}}}}}}}}}}}}}}}}}}}}}} & (5) \end{matrix}$

As for the model of the sensors 24, the model of laterally oriented sensors is provided by the expression

A _(y,meas) =ÿ _(v) +{dot over (r)} _(v) d _(xtoYA) +{umlaut over (θ)} _(v) d _(ztoRA) +r _(v) u   (6)

Accordingly, since A_(y,meas)=α_(y,body) from Equation (2), substituting the expressions for ÿ_(v), {dot over (r)}_(v), {umlaut over (θ)}_(v), and r_(v) from Equation (5) into Equation (6) yields the expression $\begin{matrix} \begin{matrix} {a_{y,{body}} = {\left\lbrack {{a_{11}{\overset{.}{y}}_{v}} + {a_{12}r_{v}} + {\frac{C_{F}}{m}\delta}} \right\rbrack +}} \\ {{\left\lbrack {{a_{21}{\overset{.}{y}}_{v}} + {a_{22}r_{v}} + {\frac{C_{F}a}{I_{z}}\delta}} \right\rbrack d_{xtoYA}} +} \\ {{\left\lbrack {{a_{41}{\overset{.}{y}}_{v}} + {a_{12}r_{v}} + {a_{43}\theta_{v}} + {a_{44}{\overset{.}{\theta}}_{v}} + {\frac{C_{F}}{m}\delta}} \right\rbrack d_{ZtoRA}} + {r_{v} \cdot u}} \\ {= {{\left\lbrack {a_{11} + {a_{21}d_{xtoYA}} + {a_{41}d_{ztoRA}}} \right\rbrack{\overset{.}{y}}_{v}} +}} \\ {{\left\lbrack {a_{12} + {a_{22}d_{xstoYA}} + {a_{42}d_{ztoRA}} + u} \right\rbrack r_{v}} +} \\ {{\left\lbrack {a_{43}d_{ztoRA}} \right\rbrack\theta_{v}} +} \\ {{\left\lbrack {a_{44}d_{ztoRA}} \right\rbrack{\overset{.}{\theta}}_{v}} +} \\ {\left\lbrack {\frac{C_{F}}{m} + {\frac{C_{F}a}{I_{z}}d_{xtoYA}} + {\frac{C_{F}}{m}d_{ztoRA}}} \right\rbrack\delta} \end{matrix} & (7) \end{matrix}$ where α_(kl) is the element in the k row and l column of the matrix A, d_(xtoYA) is the distance along the x axis from a sensor to the yaw axis, and d_(ztoRA) is the distance along the z axis from the sensor to the roll axis.

The model for vertically oriented sensors is A _(z,meas) =−g+{umlaut over (θ)} _(v) d _(yraRA)   (8) Hence, from Equations (2) and (5) $\begin{matrix} \begin{matrix} {a_{z,{body}} = {{- g} + {\left\lbrack {{a_{41}{\overset{.}{y}}_{v}} + {a_{42}r_{v}} + {a_{43}\theta_{v}} + {a_{44}{\overset{.}{\theta}}_{v}} + {\frac{C_{F}}{m}\delta}} \right\rbrack d_{ytoRA}}}} \\ {= {{\left\lfloor {a_{41}d_{ytoRA}} \right\rfloor{\overset{.}{y}}_{v}} +}} \\ {{\left\lfloor {a_{42}d_{ytoRA}} \right\rbrack r_{v}} +} \\ {{\left\lfloor {a_{43}d_{ytoRA}} \right\rbrack\theta_{v}} +} \\ {{\left\lbrack {a_{44}d_{ytoRA}} \right\rbrack{\overset{.}{\theta}}_{v}} +} \\ {{\left\lbrack {\frac{C_{F}}{m}d_{ytoRA}} \right\rbrack\delta} +} \\ {\left\lbrack {- g} \right\rbrack} \end{matrix} & (9) \end{matrix}$ where d_(ytoRA) is the distance along the y axis to the roll axis.

And for longitudinally oriented sensors, the sensor model is provided by the expression A _(x,meas) =−{dot over (r)} _(v) d _(ytoYA)   (10) such that upon employing Equations (2) and (5), Equation (10) becomes $\begin{matrix} {a_{x,{body}} = {{{- a_{21}}d_{dtoYA}\overset{.}{y}} - {a_{22}d_{dytoYA}r_{v}} - {b_{21}d_{ytoYA}\delta}}} & (11) \end{matrix}$ where d_(dytoYA) is the distance along the y axis to the yaw axis and b₂₁ is the element in the second row and first column of the matrix B.

The algorithm implemented in the estimator 26 processes the expressions from Equations (7), (9), and (11) through a filter (an estimation algorithm) to provide the estimates for the state vector {right arrow over (x)}_(v)=[{dot over (y)}_(v) r_(v) θ_(v) {dot over (θ)}_(v)]^(T).

Note that the above discussion is directed to obtaining a solution for the state vector {right arrow over (x)}_(v) in continuous time. Therefore, {right arrow over (x)}_(v), is typically discretized according to the expression {right arrow over (x)} _(v)(k+1)=A _(d) {right arrow over (x)} _(v)(k)+B _(d) {right arrow over (u)}( k)   (12) where k identifies the k^(th) time step and the matrices A and B can be discretized according to the approximations A _(d) =I _(n)+Δ_(k) ·A and B _(d)=Δ_(k) ·B where I_(n) is the nth order identity matrix, which in this case is a fourth order identity matrix, and Δ_(k) is the time step.

Although the above embodiment is directed to a sensor set with linear accelerometers, hybrid-sensor-sets are contemplated. For example, an angular rate sensor can be used in the vehicle 12 and a model of that sensor can be used in the “Kalman Filter” box 20. Specifically, for a yaw rate sensor, the model is [0 1 0 0], that is, the sensor measures yaw rate and nothing else.

Hence, in stability control, in which measuring yaw rate and roll rate/angle is useful, four accelerometers can be used for the sensors 14. Alternatively, for a hybrid system, two accelerometers and an angular rate sensor may be employed. Other examples of hybrid systems include, but are not limited to, two lateral and two vertical accelerometers; two lateral, two longitudinal, and two vertical accelerometers; and two lateral, two vertical accelerometers, and an angular rate sensor.

Other embodiments are within the scope of the claims. 

1. A system for estimating body states of a vehicle comprising: at least two sensors mounted to the vehicle, the sensors generate measured vehicle state signals corresponding to dynamics of the vehicle; a signal adjuster which transforms the measured vehicle states signals from a sensor coordinate system to a body coordinate system associated with the vehicle; and a filter which receives the transformed measured signals from the signal adjuster and processes the measured signals into body state estimates of the vehicle.
 2. The system of claim 1 wherein the filter includes a model of the vehicle dynamics and a model of the sensors, the state estimates being based on the transformed measured signals and the models of the vehicle dynamics and sensors.
 3. The system of claim 3 wherein the filter includes an estimator, an algorithm being implemented in the estimator to process the transformed measured signals and the models of the vehicle dynamics and sensors and generate the state estimates.
 4. The system of claim 1 wherein the sensors are linear accelerometers.
 5. The system of claim 1 wherein one of the sensors is an angular rate sensor.
 6. The system of claim 1 wherein the sensors include two accelerometers that measure accelerations in a first direction and two accelerometers that measure accelerations in a second direction.
 7. The system of claim 6 wherein the sensors further include two accelerometers that measure accelerations in a third direction.
 8. The system of claim 1 wherein the sensors include two accelerometers that measure lateral accelerations and one sensor measures the yaw rate of the vehicle.
 9. The system of claim 8 wherein the sensors include two accelerometers that measure the vertical accelerations of the vehicle.
 10. The system of claim 1 wherein the state estimates relate to the vehicle's lateral velocity, yaw rate, roll angle, and roll rate.
 11. The system of claim 1 wherein the signal adjuster further provides compensation for gravity biases associated with the sensors.
 12. A method for estimating body states of a vehicle comprising: generating measured vehicle state signals corresponding to dynamics of the vehicle with at least two sensors; transforming the measured vehicle states signals from a sensor coordinate system to a body coordinate system associated with the vehicle; and processing the measured signals into body state estimates of the vehicle.
 13. The method of claim 12 system of claim 1 wherein the processing includes modeling the vehicle dynamics and the sensors.
 14. The method of claim 12 wherein the generating includes measuring linear accelerations.
 15. The method of claim 12 wherein the generating includes measuring an angular rate of the vehicle.
 16. The method of claim 12 wherein the state estimates relate to the vehicle's lateral velocity, yaw rate, roll angle, and roll rate.
 17. The method of claim 12 wherein the transforming includes providing compensation for gravity biases associated with the sensors. 